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We investigate a previously overlooked bottom quark contribution to the spin-dependent cross 
section for Dark Matter(DM) scattering from the nucleon. While the mechanism is relevant to any 
supersymmetric extension of the Standard Model, for illustrative purposes we explore the conse¬ 
quences within the framework of the Minimal Supersymmetric Standard Model(MSSM). We study 
two cases, namely those where the DM is predominantly Gaugino or Higgsino. In both cases, there 
is a substantial, viable region in parameter space {m^ — < 0(100) GeV) in which the bottom 

contribution becomes important. We show that a relatively large contribution from the bottom 
quark is consistent with constraints from spin-independent DM searches, as well as some incidental 
model dependent constraints. 

PACS numbers: 12.60.Jv, 14.80.Ly 


I. INTRODUCTION 

The existence of non-baryonic Dark Matter (DM) has 
been established by many astronomical observations [I|, 
[^. Amongst the many candidates for DM, the so- 
called Weakly Interacting Massive Particles (WIMPs), 
which would have a mass in the range Oil) GeV —0{1) 
TeV, are one of the most attractive. These particles 
would only interact with Standard Model(SM) parti¬ 
cles through weak interactions (and gravity), in order 
to yield a DM relic density consistent with measurement 
= 0.1199 ± 0.0027 0. 

Direct detection of DM relies on observing the recoil 
energy after scattering from normal matter through weak 
interactions. Several DM direct detection experiments 
have claimed a possible excess, namely DAMA 0, Co- 
GeNT (sjjGRESST 0 and CDMS 0. On the other hand, 
these results are challenged by the absence of signals at 
XENONIOO 0 and LUX as well as GDMSlite 0 in 
the light DM region. The coherent, spin-independent (SI) 
interaction between a DM particle, generically labelled 
y, and a nucleus is proportional to the nucleon number. 
Because of the relatively heavy nuclei chosen for most of 
the above mentioned experiments, both the observed ex¬ 
cess and stringent exclusion limits are based on SI y — p 
scattering. 

As for spin-dependent (SD) DM detection m , in a sim¬ 
ple shell model the spin of the nucleus is that of a single, 
unpaired nucleon. As a consequence, the matrix element 
for SD y-nucleus scattering will be roughly comparable 
with that for SI y-nucleon scattering, with no enhance¬ 
ment by the nucleon number. As a result, the current 
DM direct searches place only very loose bounds on the 
SD cross section dlHil. 

In the standard calculation of SD DM-nucleon scat¬ 
tering the heavy quark contribution is usually neglected. 
That is, only the contributions from Au, Ad and As are 
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included. However, as explained in the context of the 
proton weak charge [l^, the usual decoupling of heavy 
quarks through the Appelquist-Carrazone theorem |16l | 
does not apply t o quantities influenced by the U(l) axial 
anomaly |17H2lj |. In that case, rather than being sup¬ 
pressed by inverse powers of the heavy quark mass, the 
suppression is only logarithmic. These logarithmic cor¬ 
rections were studied in considerable detail by Bass et al. 
in Refs. [H, [H, [2^, at both leading and next-to-leading 
order. As we shall explain here, there are interesting sce¬ 
narios of supersymmetry (SUSY), generally involving a 
relatively light sbottom, where the logarithmic radiative 
correction involving the &-quark that is further enhanced 
by resonant effect may make a significant contribution to 
SD DM-nucleon scattering. 

Indeed, SUSY [2^ is widely believed to provide 
the most promising explanation for new physics beyond 
SM. In SUSY models with R-parity conservation, the 
lightest supersymmetric particle(LSP) is stable and can 
become a DM candidate. On the other hand, both 
the LHG SUSY searches [2^ and naturalness argu¬ 
ments [H, [ 2 ^ suggest that only the third generation su¬ 
persymmetric quarks(squarks) can be light. In Ref. [30l| . 
it has been argued that an sbottom with a mass as light 
as ~ 0{lf)) GeV might still be consistent with current 
searches. In other models, such as the simplihed model 
framework [3l| and flavored DM models [1^ [s^ , the DM 
can only couple to the bottom quark, as motivated by 
the recent DM indirect signals . Studying the bottom 
quark contribution to the DM-nucleon SD cross section 
is crucial in models of this type. 

In this work, we focus on the Minimal Supersymmetric 
Standard Model (MSSM) with a relatively light sbottom, 
showing when and how the bottom contribution becomes 
important. When the DM is Wino, there is no coupling 
between DM and the Z-boson and only squark medi¬ 
ated processes can contribute to x-nucleon scattering. 
We investigate the parameter space where the sbottom 
contribution is comparable to, or larger than, the first 
generation squark contribution. When the DM is Hig¬ 
gsino, the first two generation squark mediated processes 
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are greatly suppressed by their small Yukawa couplings. 
However, the Higgsino can couple to the Z-boson. The 
constructive and destructive interference effects between 
Z and sbottom (b) mediated processes are discussed in 
detail for a number of variations on the structure of the 
neutralino. 

Any sbottom mediated process that contributes to the 
SD scattering cross section can also contribute to SI scat¬ 
tering. We consider the stringent LUX constraint on SI 
DM detection for light sbottom scenarios of interest. A 
relatively large SD bottom contribution can indeed be 
found, while maintaining consistency with the LUX con¬ 
straint. We also consider several model dependent con¬ 
straints from collider searches. We stress that our con¬ 
clusion has implications beyond the MSSM, which is used 
here purely for purposes of illustration. 

This paper is organised as follows. In Sec. |lTl we ex¬ 
plain the theoretical framework for the calculation of 
the SD DM-nucleon scattering cross section. Sec. IIIII 
discusses the bottom contribution for Wino and Hig¬ 
gsino DM. The corresponding SI detection and LHC con¬ 
straints on the light sbottom scenario are considered in 
Sec.IV] and Sec. El In Sec. EH we present some conclud¬ 
ing remarks. 


II. EFFECTIVE INTERACTION FOR 
SPIN-DEPENDENT DM-NUCLEON 
SCATTERING IN MSSM 


Given a general effective Lagrangian 

= dqX1^^5XQliJ.l5q , (1) 


the spin-dependent y—nucleon scattering cross section 
can be written as 


where 


p,n _ ^2| 

"30 — TV 


^SD ~ V ; J l“P,ra| 

TT rrix + mp^n 


(^p,n — 


^ ^ dq^Qp,n • 
q 


( 2 ) 

( 3 ) 


The factors A.qp^n parameterise the corresponding quark 
spin content of the nucleon: 

2s^A(jjv =< N\ipqXp"f5ij;q\N > , (4) 


where is the nucleon spin. The preferred values of the 
light quark contributions in the proton and neutron are: 

AW = = 0.84, = A^’") = -0.43, 

A(p) = = -0.02 , (5) 



FIG. 1: Processes that contribute to the Dark Matter 
spin dependent cross section for scattering from a 
nucleon. 


corresponding Feynman diagrams are given in Fig. [T] 
From those processes, we are able to calculate the coeffi¬ 
cients of the effective Lagrangian from the renormalisable 
Lagrangian below: 

C = q{aq + bqX5)xq + cq-f^x^qZ^ + dxx^XbXZp ■ (6) 

The corresponding couplings in the MSSM are written 
as [1^ 


=^%^{\z]^s^ + Z]lc^) - Y^Zj^) 

^ \/2cw 4 

Z^ 2y/2gSni ^11 y 
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zb _ 1 
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2 fjCyj 
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( 8 ) 

(9) 

( 10 ) 

( 11 ) 

( 12 ) 


We consider first the Z boson mediated amplitude i 
e non-relativistic limit: 

Msb = cd xi^i5Xj^r^^q'y’'i5q 

~ c d -\{l + 0(TO^^))x7^75X7^75g 

TTL <7 


in 


C d -^XA^lbXAiilsq 


■4'! 


where the strange quark contribution is motivated by a 
recent lattice QCD calculation [35l |. 

In the MSSM at tree level, there are two processes 
which can contribute to the effective Lagrangian. The 


so that the effective coupling dq in Eq. (HD is: 

= -£^T,q{{zUf - [Zfin (14) 
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Next, for the q mediated process we find: 

i 


■^SD = x(« - ^ 75 )g 


{Px+Pg)‘^-ml 


q{a + & 75 )x 


m? — {m^ + mqY 


— I 


m?- — {m^ + mqY 


X(a - 675)9 q{a + h'^^)x 


(a^X9 PX - b'^Xl5q 975X) 


'-x 
—z 


+ 62 


“ (™x + "1?)^ 


Xl^lhX Qliilbq ), (15) 


the mass of the 6 -quark, the q propagator in Fig. [T] can 
be effectively factored out, leaving the familiar triangle 
diagram which involves the U{1) axial anomaly. In this 
case the bottom contribution to the axial charge of the 
target proton can be taken directly from Eq. (fTTll . We 
shall take the running coupling 0:5 = 0 . 2 . As a result, 

iv) 

for the .Z-mediated process, the contribution of can 
only change the result by a factor of 


A (P) _ X\P) _ AW _ A ^ 

ffAl_ ^ 

^ A^^’) 


,{p) 


, (p) 


. (p) 


1 . 01 , 


(18) 


In this case the effective coupling in Eq. m is: 


1 0^ + 6^ 

4 m? — {rriy^ + mq)"^ 


(16) 


Note that the tree level effective coupling dq is only reli¬ 
able when rriq — is significantly larger than niq . Some 
discussions regarding the precision of the tree level ap¬ 
proximation are given in Appendix And we have also 
checked that the result calculated from Eq. dni) matches 
well with numerical tool micrOMEGAs for light flavor 
quark. 


III. LIGHT SBOTTOM CONTRIBUTION 

For most processes of physical interest the Appelquist- 
Carrazone theorem tells us that heavy quark contribu¬ 
tions are suppressed by order l/mg. However, as ex¬ 
plained in the introduction, because of the C/(l) ax¬ 
ial anomaly, the heavy quark contributions to spin- 
dependent quantities are only logarithmically suppressed. 
The particular case where this has been explored in great 
detail is the neutral weak charge of the proton. Without 
heavy quarks this is just Au — Ad — As, which has been 
used to infer values of As. However, for a precise deter¬ 
mination one must include the radiative corrections in¬ 
volving heavy quark loops which enter at order 1 / Inmg. 
For example, one finds a LO correction from the 6 -quark 
equal to [l^ : 

= - J^db(A(f) + A<f^ + AW) ^ -0.0066 . (17) 

We note that Eq. (El is second order in the strong cou¬ 
pling at the 6 mass, as is evident in the residual 5-flavor 
factor 6/23 appearing there. However, the regularisation 
of the triangle diagram leads to a logarithm in mf, in the 
numerator which has been used to cancel the logarithm 
in one factor of d&. The logarithmic radiative correction 
a[^^ is around 2 order of magnitude below the A^f ^. We 
will show later that with further enhancement from reso¬ 
nant effect the contribution from A^'^ can easily become 
dominant in the spin-dependent x-nucleon scattering. 

Provided that the difference between the sbottom mass 
and that of the DM candidate is significantly larger than 


, (p) 


which is clearly very small. 

On the other hand, the term involving A/''' can give a 
significant contribution to the spin-dependent x-nucleon 
cross section when 6 is relatively close in mass to the DM 
candidate, i.e. with resonant enhancement. For simplic¬ 
ity we study the cases where the DM particle is either 
pure Wino or pure Higgsino. The corresponding cou- 


plings are: 
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'^SD — Tvz: r~i~> y~ 


Tr^m^ + rrip 4:{m? — (m^ + rrib)^) ^ 


( 21 ) 

0.5y2(z^i)2 

( 22 ) 


where we have assumed the gauge eigenstate limit and 
only the sbottom mediated process is contributing. By 
fixing ra-^ at either 10 or 100 GeV and taking = 1 
and tan ,5 = 40, Z^ = for Wino and Higgsino DM, 
respectively, we can calculate the corresponding cross sec¬ 
tion as a function of m^. The result is shown in Fig. [2] 
From Fig. [2]we see that the sbottom can give a very large 
contribution when the mass splitting is + 100 

GeV. 


A. Comparison with the contribution from the first 
generation squarks 

First, we study the simpler case where the DM is gaug- 
ino. In this case there is no coupling between the Z bo¬ 
son and DM and only the squark mediated process can 
contribute to the SD interaction. In this subsection, we 
investigate the extent to which the sbottom should be 
lighter than first generation squark, so that they at least 
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FIG. 2: Sbottom contribution to the SD scattering 
cross section for Wino and Higgsino Dark Matter from 
the proton. 


have comparable cross sections. In the following we con¬ 
sider the sum of the contributions of all first generation 
squarks (ml, u/j, di, d/j), with their masses taken to be 
degenerate for simplicity. 

Assuming that the DM is either pure Wino or Bino, the 
ratio of the corresponding SD cross section for sbottom to 
the sum of the contributions from all the first generation 
squarks can be calculated as 




h-W 

^SD ^ _ 

qu,d-H m~ — (m^ -f 

JgD b ^ V 


,A. 


(p) 


■A 


(p) 




h-B 

qu,d — B 

^SD 


(a^ 




hi 


r? - (m^ 
. (p) 


iribY 




(p) 


- K 


(23) 


(24) 


The corresponding contours of = ag^"^ ^ are 

shown in Fig. [31 The case of Wino DM is more inter¬ 
esting than that of Bino DM because of its larger g 2 cou¬ 
pling. In this case, for 1.5 TeV first generation squarks 
and 0(100) GeV DM, an sbottom lighter than about 
200 GeV gives a larger cross section than the first gener¬ 
ation squarks. On the other hand, for Bino DM, a much 
lighter sbottom(~ 110 GeV) is required ~ too light for 
the present calculation to be reliable. 


B. Contribution coherent with that of the Z-boson 


When the DM is predominantly Higgsino, the first two 
generation squark mediated processes are greatly sup¬ 
pressed by their small Yukawa couplings. Its couplings 
to the Z boson and sbottom are dependent on the mixing 
between the two Higgsino states. 

Firstly, we briefly discuss the Higgsino mixing in the 
MSSM. In the basis (5, W, id°, id°), the neutralino mass 


mj/GeV 



m^/GeV 


FIG. 3: Contours of constant mj; which show where the 
sbottom and degenerate first generation squarks give 
the same contribution for Wino (Red solid line) and 
Bino (Blue dashed line) DM. 


matrix is given by: 

/ 


Mat = 


Ml 

0 


0 

M 2 


-cpswmz cpcwmz 
\ spswBiz —spcwinz 


—cpswmz spswBiz 
cpcwBfiz —spcwinz 

0 —/J, 

—/j, 0 

(25) 


From the mass matrix we conclude that if 


mz<.g,Mi,M2, (26) 

the four neutralino mass eigenstates Ni will be Bino 
B dominated, Wino W dominated and Higgsino (M° ± 
H°)/'/2 dominated, respectively. For example, if we also 
decouple the Bino and Wino from the mass matrix, the 
component difference between and for a given 
mass eigenstate is 


_ jjO oc 


ml 

Mig 


(27) 


For a few TeV gaugino and a few hundred GeV Higgsino, 
ANj^o_fjo ^ 0(10“^). Then, the contribution from the 
Z boson mediated process can be estimated by 


_ 12 ^ m^mp ^ 

'SD — —1-;- ) 1 / , “g 

TT rriy. + ^' 

^ q—u,d,s 




(28) 




(Af) 


TT + rrip 


, (p) 




(-i) + AWx(-i))^ (29) 


which is ~ 10 ® pb. From Fig. [31 we conclude that this 
corresponds to m-j; ~ 150 GeV for ~ 100 GeV. 
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To have a closer look at the coherent effects of Z boson 
and sbottom mediated processes, we have chosen the de¬ 
coupled Wino/Bino limit, with the Higgsino DM mixing: 

X = aUd + bHu , (30) 


where -I- = 1 and b = 1.01a, as argued previously. 

This corresponds to a cross section for the Z mediated 
process of order ~ 3 x 10“® pb. 

In this region, the b mediated process may also give a 
competitive contribution. As a result, the DM will have 
opposite sign coherent effects for the proton and neutron: 


X-iP,n) _ 
^SD “ 


12^ m^TOp 


TT rriy- + rup 
+ +T3sA<f’^^)- 


8m^ 

mt)^) » 


4(to? — (m 


(31) 


This makes the detailed consequences for SD DM scat¬ 
tering from real nuclei potentially very complex. 

We show the importance of the 6-quark contribution 
through its coherent effects between Z mediated and 
sbottom mediated processes in Fig. 01 There we have 



0 100 200 300 400 500 600 700 800 900 1000 


mj/GeV 
= 100 GeV 



mj/GeV 


FIG. 4: Constructive and destructive coherence effects 
between Z mediated and b mediated processes for 
= 10 GeV (Upper) and = 100 GeV (Lower), 
respectively. The vertical arrowed lines at mi = 30 GeV 
and mi = 120 GeV indicate the value below which the 
calculation should not be considered reliable. 

marked out the compressed spectrum region [Aimi — 
™x) ^ 20 GeV) where our tree level calculation cannot 


be considered reliable. According to the Eq. with 
assumed Z'^ = I.OIZ^^, the first term in the parenthesis 
is positive for proton and negative for neutron, while the 
second term is always negative because < 0. 

As a result, the interference terms for proton and neu¬ 
tron are constructive and destructive, respectively. We 
find that for the DM mass around 0(10 — 100) GeV an 
sbottom with mass < 300 GeV can make a non-negligible 
contribution. In some specific regions, the corresponding 
cross section for the Z-mediated process may even be en¬ 
hanced or reduced by several orders of magnitude. 


IV. SPIN-INDEPENDENT DM DETECTION 
CONSTRAINT FROM LUX 

The same process shown in the right panel of Fig. [l] 
which can give rise to an enhancement of the spin- 
dependent scattering cross section, can also contribute to 
spin-independent scattering. As a result, the very strin¬ 
gent spin-independent DM search bound from LUX 0 
may already exclude some of the parameter region found 
to be of interest here. 

We start with the following effective Lagrangian [3^ 

0 : 

„(i) 

^Si = ^{fqmqXXm+^Xid^YxOl^ 

q X 

( 2 ) 

+ ^-\x{^^>^){^^nxO%) + (32) 

where x is the DM field, its mass and the twist-2 
operator: 

~ ~ ■ (33) 

The corresponding spin-independent scattering cross sec¬ 
tion of DM with a proton can be written as 

^sV = (34) 

TT 

where /i = m-^m^/{m^ -|-mjv) and 

~ = 4(9(2)+ 9(2))(4^^+4^^) 

q—u,d,s q—u,d,s,c,b 

^ / 9 /U+ ^(9(2)+ g(2))(4^^+4^^) 

q—u,d,s q—u,d,s,c,b 

+ ^ X! fTcfQ (36) 

Q=c,b,t 

The light quark parameters /t, are defined by 

=< N\mqqq\N > , (37) 
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and /tg = 1 — J2q=u d s fTq- Recent lattice simulations 
give Ulil: 

fP = 0.023, /^ = 0.033, /f = 0.026. (38) 


The second moments of the parton distribution func- 
tions(PDFs) can be used to evaluate the matrix element 


- ^m%gij_^){q{2) + q{2)) = ttin < N{p)\Ol_jN(p) > , 

(39) 

which from the CTEQ PDFs [i^ yields 

6(2) = 0.012, 6(2) = 0.012 , (40) 


at the Z boson mass scale. 

Using a similar technique to that used in calculating 
the SD effective coefficient, dq, above, we can find the cor¬ 
responding effective coefficient for the spin-independent 
case. Based on the renormalizable Lagrangian Eq. ®, 
we have 


f,= 


- {rux + rUqYY 8 


a, - 


mq{m^—rUqY) 4 


g^q^ + gf^ = (_2 


ruy 


« (m? - (mx+m5)2)2 2 


(41) 

(42) 


As a result the twist-2 operator, gives a much larger 
contribution than fq in most cases. For Higgsino DM, 
where aq ^ bq, the second term of fq can easily become 
dominant. However in this case it is negative, so a can¬ 
cellation between fq and Pq may happen in some of the 
parameter regions. 

We first consider the pure Wino DM case, with only 
6z, mediated scattering. From Eq. (l36)) . we have 


+ b{2)){gl--^) + (43) 

where 


rbL-w ^ 1 

32 (m? - m2 )2 

bL-W _ g^T^X 1 

^ 8 (m~ — m2)2 


(44) 

(45) 


The tree level calculation for the spin-independent and 
spin-dependent cross sections is shown in Fig. [5] We 
conclude from the figure that > 50 GeV is 

required to evade the spin-independent bound from LUX 
for Wino DM. It has to be noted that the pole at m-^ = 
rrib -b m^, for SI tree level results, will not show up when 
the full NLO effects are taken into account . We have 
checked that our results fit the numerical results from 
micrOMEGAs [11,14^ quite well, outside the pole region. 


-Q 

Q. 


b 
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mj/GeV 


FIG. 5: SD and SI scattering cross section for Wino 
DM from the proton, with m^ = 10 GeV and = 100 
GeV, respectively. 


Next, we discuss the more interesting case where the 
DM is predominantly Higgsino. As discussed above, in 
this case the relatively large spin dependent cross sec¬ 
tion from the sbottom mediated process can interfere co¬ 
herently with the Z mediated process, leading to very 
different SD scattering rates for protons and neutrons. 

The SI DM-proton effective coupling is 


= mp( J(6(2) + h{2)){gl--^) + ^fTaft~^) 

where 

fbi-H ^_Wx_ 0.5Ya^(Z|/)^ 

^ (m? - {m^ + mb)^)"^ 8 

_^_ Y,^{ZffZlZ^ 

mb{m? - {m^ + mb)'^) 4 

^ [m? - {m^ + irib)^)^ 2 


(46) 


(47) 

(48) 


The corresponding tree level calculation for the spin- 
dependent and spin-independent cross sections is shown 
in Fig. [6] From that figure we see that a small compo¬ 
nent of left-handed sbottom is favoured in order to evade 
the LUX bound. When the left-handed sbottom com¬ 
ponent is relatively large, a large SD cross section may 
also be consistent with the LUX experiment if there is a 
cancellation in crsi- 


V. MODEL DEPENDENT CONSTRAINTS AND 
A GENERAL ARGUMENT 

We have presented a representative study of the poten¬ 
tial importance of the bottom quark contribution to DM 
spin-dependent detection within the framework of the 
MSSM. This particular contribution has hitherto been 
overlooked. However, in a realistic model such as MSSM, 
there are many other incidental constraints. We will 
briefly outline how these may be evaded, while keeping 
our discussion as general as possible in this section. 
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0 50 100 150 200 250 300 



0 50 100 150 200 250 300 


mi! GeV 

FIG. 6: SD and SI scattering cross section for different 
mixing of Higgsino DM from the proton. Upper: 
rriy. = 10 GeV, Z~^ = 0.01. Lower: = 10 GeV, 

= 0.04. 

0 

LEP placed a very stringent bound on the chargino 
mass > 92.4(91.9) GeV) [s^. Because for 

either Wino and Higgsino DM there is a charged partner 
(chargino), which has very similar mass with the DM, we 
cannot have Wino and Higgsino DM of < 90 GeV in 
a typical MSSM framework. 

As for rriy^ > 100 GeV, on the other hand, it will be 
constrained by LHG sbottom searches [HI, and mono¬ 
jet search [^, since we usually need a relatively light 
sbottom to enhance the bottom quark spin dependent 
contribution. The corresponding LHG exclusion bounds 
and spin-dependent x~ P scattering cross section in the 
mjj — plane are shown in Fig. [T] To generate this 
figure we have used Eq. (I^Tl) and (0^ . where only the 
sbottom mediated process is considered. The contours 
of show the condition when the sbottom mediated 
contribution is half the size of Z mediated process for 
Wino DM and a typical Higgsino DM candidate with 
Zf^ = I.QIZ^. This figure suggests that a large por¬ 
tion of parameter space is excluded by the LHG sbottom 
searches. 

However, there are several ways to avoid these con¬ 
straints: 

• For rriy^ < 100 GeV, we can work in a more general 
framework, where the dark matter does not have 
any charged partners. Its couplings to the Z boson 
and b may be of the same order; e.g. the simpli¬ 



FIG. 7: LHG exclusion bound on sbottom mass versus 
DM mass. The grey, green and cyan shaded region 
correspond to the exclusion limits given by LHG 
searches for sbottom at 8 TeV 1^. 1 3 TeV [5l| and 
searches for mono-jet at 8 TeV [^, respectively. 

fied model framework [Hl | or flavored dark matter 
models [H, js^. 

• For 100 GeV, if the charged Higgsino decays 

into DM and a relatively long lived particle, with 
lifetime of 0(10—100) cm, similar to Ref. [sj. As a 
result, the reconstructed track will not point to the 
interaction point and would therefore be unlikely to 
be considered a ’’good” track. In this case, the LEP 
constraints on charginos can be evaded. The light 
sbottom constraint can also be evaded by tuning 
appropriate mixing ~ see e.g. Refs. [13, 

• For > 100 GeV, if the sbottom is decayed in 

more complicate modes other than b —> the 

corresponding LHG bound on sbottom mass can 
be loosened. 

• We can also work with heavier DM, e.g. m-^ = 300 
GeV, for example, as shown in Fig. [8] In this case, 
^6 ^ 350 GeV is consistent with the LHG searches, 
while the sbottom mediated process can give a sig¬ 
nificant contribution to spin-dependent scattering. 


VI. CONCLUSION 

In this work, we have demonstrated the potential im¬ 
portance of the bottom quark contribution to the DM 
spin-dependent cross section due to the axial anomaly 
and resonant enhancement, which has hitherto been over¬ 
looked. Even though our calculation was carried out 
within the framework of the MSSM, the general conclu¬ 
sion will be relevant to any models with similar particle 
content, since the only relevant ingredients are y, b and 
the Z-boson. 

In the MSSM, we calculated the bottom quark contri¬ 
bution to spin-dependent x ~ P scattering. Firstly, we 
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FIG. 8: Constructive and destructive interference 
effects between Z mediated and b mediated processes 
for Higgsino DM with = 300 GeV. 


considered Gaugino DM, where there is no coupling be¬ 
tween the Z-boson and DM. Assuming = 100 GeV 
and degenerate first generation squarks with a mass of 
1.5 TeV, we found that an sbottom of mass < 200 
GeV can give rise to a larger spin-dependent cross sec¬ 
tion for Wino DM. By contrast, for Bino DM a much 
lighter sbottom mass (mj, ^110 GeV) is required to give 
a competitive cross section. As for Higgsino DM, the first 
generation squark contributions are suppressed by their 
small Yukawa couplings. However, the Z-boson mediated 
process does contribute. For a given Higgsino mixing of 
DM the sbottom mediated process may interfere either 
constructively or destructively with Z-boson mediated 
processes, with different signs for protons and neutrons. 
For a typical mixing of Higgsino DM with mass around 
0(10— 100) GeV, we find that an sbottom of mass below 
300 GeV can have non-negligible effects. 

The squark mediated process that gives rise to an in¬ 
crease in spin-dependent DM scattering can also con¬ 
tribute to the spin-independent cross section. Our cal¬ 
culation shows that A(mj, — rn^) > 50 GeV is required 
to evade the LUX constraint for Wino DM, while for 
Higgsino DM, either a small component of left-handed 
sbottom or a large cancellation in ctsi is needed. Some in¬ 
cidental model dependent constraints from LEP and the 
LHC are considered as well. Those constraints, however, 
can be evaded in more general theoretical frameworks. 

As pointed out earlier, our tree level results may break 
down as the sbottom and DM masses become degenerate 
(A(mj; — m^) < 20 GeV). We leave the higher order 
calculation for this small region for future work. Finally, 
while the calculations for the top quark case will be more 
complicated because there is no clear separation of mass 
scales for interesting ranges of DM mass, there is a clear 
need to investigate the role of the axial anomaly for that 
case too. 
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Appendix A: Precision of tree level approximation 



The heavy quark contributions to the axial charge 
start at two loop level through the process shown in 
Fig. m A detailed calculation of this diagram is given 
in Ref. [H, [2^. In this study, all we need to know is 
the rriq dependence of the amplitude. Then we can de¬ 
rive the range of mg for which the tree level effective 
coupling Eq. [iniis justified. 

The mg dependence only exists in the triangle loop 
that is marked by the grey shaded ellipse in Fig. |9l The 
vertex amplitude is 


-pab _ 

^ fJ-aP ~ 



rp r 7 a»75 

( 2^)4 


' — m I 






(Al) 


where m,q is the quark mass and momentum in the tri¬ 
angle, is the dark matter four momentum and mg is 
the squark mass. After introducing Feynman parameter 
X and the substitution: 


P = , (A2) 

A = x^p^ — xp^ + m'^{l — x) + xmg , (A3) 

the vertex amplitude Eq. ED can be simplified to 

/■i r iPi^ 

oc 3! dx {-2x{p^)a J ~ 

f dH P ,,22 2^ /■ dH 1 

J (2^)4 (^2 _ A)4 Px "4 ) y (27r)4 (;2 _ a)4 ' 

(A4) 
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FIG. 10: Upper: The uia dependence of the ratio for 
several given DM masses. Lower: In the shaded region, 
the tree level approximation matches the loop level 
results within required precision. 


On the other hand, the rriq dependence of the tree level 
effective coupling can be factored out as 


oc 


[m^ 


1 

)2 — rn'~ 


(A6) 


So we can define the ratio 


Ratio oc r'°°p/r‘™ 

= ((m^ - mir - ml) ■ . (A7) 


Taking the non-relativistic limit for the DM momen¬ 
tum, i.e. = (m^, 0,0, 0), m = mf, and mq = + 

we solve the ratio numerically. The results are shown in 
the upper panel of Fig. [TO] From the figure we find that 
the ratio tends to a constant in the heavy squark region 
for given m^^ which means the tree level description is 
accurate. While in the region of small mass splitting, the 
tree level results deviate from the full loop calculation 
considerably by a amount. 

The range of tua that permits the tree level approx¬ 
imation in required precision P can be solved by using 
the inequality 


Ratio(mA)/Ratio(mA = 500) > (1 — P) (A8) 


Finally, after integrating out the P, we will get a simple 
mq and p^^ dependence of F“^^: 





xpf^^x^p^. — m?) 
A2 


) (A5) 


at each given DM mass. In the lower panel of Fig. [TOl 
we show the tua region in which the tree level approx¬ 
imation matches the loop level result within 20% and 
50% precision, respectively. For example, when m^ ^ 10 
GeV, ruA > 20 GeV is sufficient to guarantee that the 
tree level approximation is accurate within 20%. 
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